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The structural characteristics of the perovskite-based ferroelectric Pb(Zni/3Nb2/3)i-xTix03 
at the morphotropic phase boundary (MPB) region (x~ 0.09) have been analyzed. The analysis 
is based on the symmetry adapted free energy functions under the assumption that the total 
polarization and the unit cell volume are conserved during the transformations between various 
morphotropic phases. Overall features of the relationships between the observed lattice constants 
at various conditions have been consistently explained. The origin of the anomalous physical 
properties at MPB is discussed. 
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§1. Introduction 

The physical properties in solids at the phase bound- 
aries between the two competitive stable phases tend 
to exhibit anomalous characteristics. A typical example 
is seen in the extraordinary transport phenomena such 
as high-Tc superconductivity and colossal magnetoresis- 
tance at metal-insulator phase boundary of the transi- 
tion metal oxides. The high piezoelectric effect at the 
morphotropic phase boundary (MPB) observed in sev- 
eral perovskite-based ferroelectric materials may be con- 
sidered to provide another interesting example of such 
phenomena at the structural phase boundary. 

As is discussed in Part I^^ of the present paper, the 
high piezoelectricity has been investigated from both ex- 
perimental as well as theoretical aspects. Experimen- 
tally, Noheda et al?^^^ discovered a previously unknown 
monoclinic phase belonging to space group Cm (Ma- 
phase) in the boundary region between rhombohedral 
(R) and tetragonal phases in PbiZri-xTixOa (abbrevi- 
ated as PZTx). This is immediately followed by the the- 
oretical work by Bellaiche et al..^^ Based on the first 
principles calculations they showed that there should in 
fact exist Ma-phase within a narrow concentration range 
intervening between the high symmetry R- and T- phase 
in PZTx system in complete agreement with the experi- 
mental results. It is also pointed out that the high piezo- 
electricity is associated with the stabilization of the mon- 
oclinic phase. 

Subsequently, X-ray measurements by Noheda et alJ^ 
on another type of material Pb(Zni/3Nb2/3)i-xTix03 ( 
abbreviated as PZN-xPT ) also revealed the existence of 



P. 




Fig. 1. The trajectory followed by Ps in the 3-dimensional po- 
larization vector space. R, T and O represent the symmetries of 
the structure as Ps is oriented to [111], [001], and [101] respec- 
tively. The symbols A and B indicate the process between R— >T 
and T— >R respectively. 



similar monoclinic phase at the MPB region, except that 
the space group in this case is Pm (Mc-phase) instead of 
Cm. 

On the other hand, Ishibashi and Iwata^"^°) developed 
a phenomenological theory on the physical properties at 
morphotropic phase boundary in perovskite solid solu- 
tions. He pointed out that depending on the coefficients 
of fouth order terms of the free energy expansion, the 
system stabilizes the orthorhombic phase intervening be- 
tween the tetragonal and the rhombohedral phases. He 
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also discussed the dielectric*' as well as elastic^") re- 
sponses of the system at MPB and pointed out that the 
anomalous behavior in the susceptibiUty tensor compo- 
nents is in fact due to large electromechanical coupling 
constants of these materials. Fu and Cohen^^^ presented 
the key concept to understand the physics of high piezo- 
electricity: that is, the rotation of the spontaneous po- 
larization during the structural transitions under electric 
field. Through first principles calculations they showed 
that there is a low energy path along the lines combin- 
ing the symmetric directions of [111], [001] and [101] as 
shown in Fig.l. Upon application of the external field, 
the polarization would easily rotate along the particu- 
lar path, which results in a large piezoelectric response. 
The stabilization of the monoclinic phase are considered 
to manifest this view point since monoclinicity is induced 
during the polarization vector continuously rotates along 
the specific path. 

As is discussed by Noheda et al.J^ in the case of PZN- 
8PT, the Mc-phase appears on the path between P \\ 
[001] and P \\ [101]. This suggests that the high sym- 
metry stable phase for composition around 8PT would 
be tetragonal and orthorhombic (O) rather than rhom- 
bohedral. In fact, depending on the prehistory of the 
sample preparation such as heat treatment, poling, etc., 
the orthorhombic phase is observed to become stabilized 
in PZN-9PT powder sample.^^) 

Recently, we carried out a neutron diffraction study 
using single crystals of PZN-9PT, with focus on the ex- 
istence of Mc-phase. The results are summarized in Part 
I of the present paper. In the present study, we develop a 
theoretical treatment based on a phenomenological argu- 
ments in order to analyze the experimental results pre- 
sented in Part I. Since the treatment is phenomenolog- 
ical, unlike the treatment based on the first principles 
calculations, the conclusions depend on a few disposable 
parameters included in the free energy expression. Nev- 
ertheless, the results seem to provide some important 
insights into the specific physics at MPB regions. 

§2. Analysis of Experimental Data 

Using BaTiOa as a model system, Fu and Cohen^^^ 
developed a detailed theoretical treatment and dis- 
cussed that the morphotropic phase transition scheme in 
BaTiOs is described by the 'rotation' of the spontaneous 
polarization vector, Pg , within the {P^, Py, Pz) space 
moving along [111]^[001]^[101] directions as shown in 
Fig. 1. This process realizes the phases with rhombo- 
hedral, tetragonal and orthorhombic symmetry. Later, 
Vanderbilt and Cohen^^^ pointed out that when the en- 
ergy terms up to the eighth order with respect to po- 
larization are included in the free energy expansion, the 
monoclinic phases with the polarization along [^ ^ Q- 
(Ma-phase) and [^ C]-(Mc-phase) directions can be sta- 
bilized. 

Following these previous arguments, we also assume 
that the polarization moves on the surface of the sphere 
with |Ps| conserved during the process of the polymor- 
phic phase transitions in the present system. For later 
convenience, let us define the parameters to specify the 
orientation of polarization as it moves along the trajec- 



tories A and B in the P-space (See Fig.l.) as follows: 
process A: 



P=p2 



process B: 



P = 



p2 



(2.1) 



(2.2) 



In the process A, p changes in the range of i < p < 1 
as P rotates from [111] (p = i) to [001] (p=l), while in 
the process B, p' changes in the range of < p' < | as 
P rotates from [001] (p'^O) to [101] (p' = i). 

In addition, we further assume that the volume of 
the unit cell is also conserved throughout the sequen- 
tial phase transitions. As is typically seen in the case of 
BaTiOa, this assumption seems to be satisfied well. 

Since we are particularly interested in the lattice dis- 
tortions, we give the expression of the 'elastic' free energy 
leaving the strain components as the independent vari- 
ables rather than projecting on the polarization space. 
In order to make the treatment physically transparent, 
we construct the basis functions of the six dimensional 
direct product space of {P^,P^,P^,PyPz,PzPx,PxPy) 
and {cxx, Syy, P-zz, P-yzT Szx, Sxy) associatcd with the irre- 
ducible representations of the point group m3m as fol- 
lows. 



i;, = —{P^ + P^ + P^)eA^, 



^, = ^(2P^-P^-P^) 



and 



V6 

V'4 = PyPz 

i'5 = PzPx ) e T2g 
tp6 = PxPy 



(2.3) 



62 — ^(26x2 Cxx Gyy) 

63 = 775 (Sxx ~ ^yy) 



e Ea 



V2 

64 = 6y, 

65 = Szx ^ e T2g 

66 ~ ^xy 



(2.4) 



Then, we can set up the simplest symmetry adapted elas- 
tic free energy including only the quadratic (or bilinear) 
terms as follows: 



where 



TP _ p('4ig) , Tj{Eg) 771(729) 



(2.5) 

(2.6) 



Fef'^ = ^(4 + 4) - fl'(e2V'2 + 63^3), (2.7) 
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(b) Process B 



\c"{el + el + el) 



-g"{e4'4'4 + e5tlj5 + eetpe). (2.8) 



The second terms of these equations give the electrostric- 
tive energies in the cubic perovskite system. 

In these expressions, the independent thermodynami- 
cal variables are e^s {v = 1, • • • 6) while tjy^s are consid- 
ered to be the parameters to define the orientation of 
the polarization which are explicitly given by p and p' 
in eqs. (2.1) and (2.2). Among the three terms in eq. 
(2.5), F^f^^ is irrelevant to determine the stability of 
the system during the rotation of the polarization vec- 
tor, since it is expressed in terms of ei and ipi which are 
both assumed to be conserved. 

(E \ 

Moreover, F^^ " is defined in the 2-dimensional space 
spanned by the basis of Eg representation (e2 and 63), 
while F^f^"^ in the 3-dimensional space spanned by the 
basis of T2g representation (64, 65, and ee)- Therefore, in 
order to find the stable lattice distortions, we simply min- 
imize F^f"^ and F^f^''' independently. 

(i)£Jg -space 

By the standard minimization procedure, it is easily 
shown that the trajectory of the minimum energy €^{p) 
is given by: 

(a) Process A 



Hp) 



0. 



1), 



(2.9) 



9 (I . 



(2.10) 



where £2 and ea are defined by 120° rotation of the (e2,e3) 
coordinate system. 

The trajectory followed by the representative point of 
the system in the (£2,63) space is given in Fig. 2(a). In the 
figure, R, T and O represents the rhombohedral, tetrag- 
onal and orthorhombic phases respectively. 

(ii)r2g-space 

Similarly, we have the trajectory of the energy minimum 
path: 

(a) Process A 



el{p) = el{p) 



(b) Process B 



^1 = 4. = ^. 

4{p')^'-^p'i^-p')- 



(2.11) 



(2.12) 



The trajectories followed by the representative point 
of the system in the (64, es, eg) space is given in Fig. 2(b). 

Eqs. (2.9) to (2.12), together with the condition of the 
volume conservation (e°^ + -|- e°_j, = 0), give the set 
of lattice parameters of the stable state in terms of p or 
p' as shown in Table I. Notice among total of 30 lattice 
constants, there are only two disposable parameters, A 
and 5, besides p and p' . We utilized the values of c\ (the 
unique axis in T-phase) observed in 9PT5 and /3 (mon- 
oclinic angle in 0-phase) in 9PT4 as tabulated in Part I 
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Table I. The lattice distortions of the perovskite system in the polymorphic phases realized when Ps rotates through the process A 
and B. The disposable parameters are defined by: A = ^t,S = ^77. 
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(sec Tabic I) to obtain the value of A and 5 respectively. 

Fig. 3 shows the sequence of the calculated lattice con- 
stants to be taken by the system as p and p' is varied 
to stabilize (R)-^(Ma)^T^Mc^O, including the ficti- 
tious stable phases R and Ma which are not realized in 
the PZN-9PT system, but realized in the related mate- 
rials. The observed values at lower temperatures given 
in Part I are indicated by the solid circles. We can not 
calculate the lattice parameters as a function of tempera- 
ture for direct comparison with the experimental results, 
since the temperature dependences of p{T) and p'{T) are 
outside of the framework of the present treatment. 



O 



Mc 



(Ma) 



(R) 




Fig. 3. The calculated variation of the lattice constants in PZN- 
9PT as the parameters p and p' is varied. The solid circles are the 
experimental values obtained by the present work (Part I). The 
figure includes the fictitious phases (R and Ma) which are not 
realized in 9PT, but realized in the closely related systems. The 
open squares indicate the data used for the numerical analysis. 
(See the Text.) 



§3. Conclusions and Discussions 

In conclusion, based on the symmetry adapted free 
energy function, we analyze the structural characteris- 
tics of PZN-9PT imder the assumption that the total 
polarization and the unit cell volume are conserved dur- 
ing the transformations between various morphotropic 
phases. Overall features of the relationship between the 
observed lattice constants in various phases have been 
consistently explained. 

So far, we have discussed the lattice distortions which 
are directly observed by the diffraction studies described 
in Part I. The energetic considerations concerning the 
stability of the phases seem to give more important in- 
sight into the physical properties of the system at MPB 
as discussed in the following. 

Using eqs. (2.9) through (2.12), the free energy values 
of the rhombohedral (Fe((R-))> the tetragonal- (i^ei(T)), 
and the orthorhombic-(Fe;(0)) phases are given by, 



FeliO) 



3c' ' 

12c' 



g 

8c" ' 



(3.1) 



which means that (Fe;(0)) is expressed by a fixed 
weighted mean of (Fe;(R)) and {Fei{T)) as, 

F,i{0) = \f,i{T) + ^F,i{R), (3.2) 

irrespective of the parameter values of c's and g's. As- 
suming, for simplicity, linear dependence of the param- 
eters on the x- value (the concentration of PbTiOs), the 
relative energies of the three phases change with respect 
to X as schematically shown in Fig. 4. We see that the sta- 
ble phases are generally given by either T- or R- phases 
while 0-phase barely becomes one of the three degener- 
ated ground states at the single critical concentration, 

^C5 

where, 



r,'2 



9"' _ 



(3.3) 



3c' 6c" 

is satisfied. Notice the critical concentration x = x,_ 
nothing but the MPB point. (See Fig.4.) 

More generally, the free energy Fei{p) and Fei{p') with 



IS 



Symmetry of high-piezoelectric Pb based complex perovskites at the morphotropic 



5 




Fig. 4. Schematic diagram of the relative phase stability as x- 
value (the concentration of PbTiOs) is changed. The orthorhom- 

bic phase is stable only at the single critical concentration, Xc, 
which just corresponds to MPB point. 



an arbitrary polarization direction is given by, 

„'2 



-|^(1-P)(1 + 3P), 



(3.4) 



Fel{p') = - 



12c' Ac 



7(1 -V)^ 



(3.5) 



It is easily seen that at MPB where eq. (3.3) is satisfied, 

F,i{p) = Fei{p') = A, (3.6) 

irrespective of the values of p and p' which means that at 
MPB, the energy is completely degenerated throughout 
the postulated path of polarization rotation. Only by the 
introduction of higher order energy terms which are ne- 
glected in the present treatment, the degeneracy will be 
lifted so that either of 0-, Ma- or Mc- phases may become 
stabilized in the MPB region as shown by Vanderbilt and 
Cohen. This situation may explain the experimental 
observations that depending on subtle differences of the 
prehistory of sample preparation, the symmetry of the 
stable state becomes different. 

Furthermore, neglecting the higher order terms, we 
see that at x = Xc the system could move over between 
R<^=^Ma<J4>T<^Mc<;4>0 under external forces such as elec- 
tric field and stress without any cost of elastic energy. 
Therefore, at MPB some of the susceptibility tensor com- 
ponents of Xij (dielectric constants), dijk (piezoelectric 
constants) and the inverse of Cijki (elastic constants) 
should critically blow up, although in the real system, 
the higher order effects would tend to prevent such a 
strong anomaly. Actually, Ishibashi and Iwata^"^°) have 
given the explicit expressions for these susceptibility ten- 
sors, some of which tend to show the critical behavior. 
The above simple energy considerations seem to reveal 
the essential origin of the anomalous physical properties 



at MPB. 

As a further application of this view point, we take 
notice on the elastic anomaly exhibited by hexagonal 
BaTiOa.") It is known that hexagonal BaTiOs (h-BT) 
undergoes two phase transitions successively as hexago- 
nal orthorhombic (space group C222i) monoclinic 
(space group P2i).^'''' The relevant order parameter 
is the doubly degenerated phonon modes, Qi and Q2. 
Ishibashi^^' -^^^ pointed out that the transition between 
O- and M- phases is caused by the rotation of the polar- 
ization vector of the mode in the 2-dimensional (Qi,Q2)- 
space, which is very similar to the feature presented in 
this paper if we replace P by Q. However, it should be 
noticed that the anomaly at the 0-M transition region 
in h-BT is manifested in the elastic compliances rather 
than piezoelectricity because the Q-modes are optically 
inactive. In this context, h-BT may be considered to pro- 
vide another example belonging to the different category 
concerning MPB adventure. 
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